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Abstract 

A gamma-ray burst (GRB) is a strong and fast gamma-ray emission from the explosion of stellar systems (massive stars 
or coalescing binary compact stellar remnants), happening at any possible redshift, and detected by space missions. 
Although GRBs are the most energetic events after the Big Bang, systematic search (started after the first localization 
in 1997) led to only 374 spectroscopic redshift measurements. For less than half, the host galaxy is detected and studied 
in some detail. Despite the small number of known hosts, their impact on our understanding of galaxy formation and 
evolution is immense. These galaxies offer the opportunity to explore regions which are observationally hostile, due 
to the presence of gas and dust, or the large distances reached. The typical long-duration GRB host galaxy at low 
redshift is small, star-forming and metal poor, whereas, at intermediate redshift, many hosts are massive, dusty and 
chemically evolved. Going even farther in the past of the Universe, at z > 5, long-GRB hosts have never been identified, 
even with the deepest NIR space observations, meaning that these galaxies are very small (stellar mass < 10^ Mq). 
We considered the possibility that some high-z GRBs occurred in primordial globular clusters, systems that evolved 
drastically since the beginning, but would have back then the characteristics necessary to host a GRB. At that time, the 
fraction of stellar mass contained in proto globular clusters might have been orders of magnitude higher than today. Plus, 
these objects contained in the past many massive fast rotating binary systems, which are also regarded as a favorable 
situation for GRBs. The common factor for all long GRBs at any redshift is the stellar progenitor: it is a very massive 
rare/short-lived star, present in young regions, whose redshift evolution is closely related to the star-formation history 
of the Universe. Therefore, it is possible that GRB hosts, from the early Universe until today, do not belong to only one 
galaxy population. 
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1. Introduction 

The discovery of a galaxy hosting a gamma-ray burst 
(GRB) was achieved for the first time in February 1997, 
with the identification of an afterglow, the event GRB 970228 
(Gosta et all997). However, its redshift, z = 0.695, was 
measured later from the detection of emission lines in the 
host through Keck spectroscopy (Djorgovski, et al. 1999). 
The very first gamma-ray burst redshift was measured in 
May 1997, for GRB 970508. The absorption lines seen in 
the Keck spectrum of the optical afterglow gave z = 0.835 
(Metzger et al. 1997). Today, precise localization of after¬ 
glows (2 arcsecs or better) is routinely performed, at the 
level of several events per week, mainly thanks to the data 
collected with the most successful dedicated space mis¬ 
sion, the NASA satellite Swift, launched at the end of 2004 
(Gehrels et al.2004). Since 1997, and as of the first half 
of 2015, gamma-ray instruments identified a total of more 
than 1,400 GRBs. Although ^ 90% of Swift GRBs are 
localized, thanks to the X-ray instrument XRT, the after¬ 
glow localization for the whole population since 1997 was 
possible for only about half them. These are mainly long 
GRBs (more than 90%), those for which the gamma-ray 
emission is longer than a couple of seconds, and associated 


with the death of a massive star (mass of the progenitor 
M > 30 Mq). The hosts are detected mainly for the long 
GRBs. 

Fruchter et al. (2006) investigated the location where 
the GRB takes place and found that those at z < 1.2 
prefer the most active regions in the galaxy, more than 
what done by core-collapse supernovae (GG-SNe). Follow¬ 
ing this early result, Kelly et al. (2008) found that GRB 
environments are more similar to those of SN type Ic rather 
than SN type II. This is consistent with the fact that spec¬ 
troscopically confirmed SNe found a couple of weeks after 
the GRB are type Ic (see Hjorth & Bloom 2012; for a 
review). SN Ic’s tend to be more luminous than the typi¬ 
cal GG-SNe. Studies of the SN-GRB connection is limited 
basically to z < 1, where SN can spectroscopically be iden¬ 
tified. 

At higher redshift, it has been finally found that the 
role of dust is important. The TOUGH sample includes a 
complete investigation of 69 GRB hosts (median redshift 
z = 2.14 ± 0.18). Hjorth et al. (20112) indicate that opti¬ 
cally dark GRBs tend to occur in more massive and red¬ 
der galaxies with respect to optically detected GRBs. The 
presence of large dust content was found in dark GRBs 
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Figure 1: Number of GRBs per redshift bin (as of May 2015). All 374 redshifts are spectroscopically determined, either from the afterglow 
or from the host galaxy. The blue curve and yellow shaded area represent the star-formation rate density (SFRD) of the Universe (from Wei 
et al. 2014), scaled to the GRB histogram for 1 + z < 2. The SFRD below 1 + 2 : ~ 5 is the one determined by Hopkins & Beacom (2006) and 
Li (2008) from an observational compilation of UV galaxies. The 1 + 2 ; > 5 SFRD (yellow shaded area) is vaguely constrained by the GRB 
detection rate (Ghary et al. 2007; Yiiksel et al. 2008; Wang & Dai 2009). 



though radio and deep observations. The large sample 
SHOALS studied by Perley et al. (2015a,b) with multi¬ 
band optical-IR observations of 119 hosts revealed a rela¬ 
tively large fraction (20%) of dust obscured galaxies, which 
are also massive systems. Radio observations (Michalowski 
et al. 2012; Perley & Perley 2013) have revealed that some 
have properties similar to sub-millimeter galaxies, but the 
majority does not have hidden high dust obscured star- 
formation rate (SFR). 

Despite the numerous surveys which used with different 
means and selection criteria, we are still not sure how well 
GRB hosts represent the whole galaxy population. For 
instance, contradictory results are found on the galaxy lu¬ 
minosity function. Using the THOUGH sample, Schulze 
et al. (2015) compared the GRB-host luminosity function 
to the one of Lyman break galaxies (LBGs) and concluded 
that GRBs select metal poor galaxies. This was not con¬ 
firmed by the high-redshift 3 < z < 5 sample studied by 
Greiner et al. (2015). 

Since the beginning, it was pursued the idea that GRBs 
prefer low-metallicity environments (e.g., Graham & Fruchter 
2013; Vergani et al 2014; and references therein). From the 
theoretical and modelling point of view, this conclusion 
was reached at low (Niino et al. 2011; Boissier et al. 2013; 
Vergani et al. 2014) and high redshift (Ghisari et al. 2010; 
Trent! et al. 2015). However, see the work by Gampisi et 
al. (2011) and Elliott et al. (2012) for a different conclusion. 

This shows that we are still far from a full comprehen¬ 
sion and interpretation of the galaxy population hosting 
GRBs. Gontrary to what is commonly believed, our main 
limitation is not the observational bias, but the small num¬ 
ber statistics. Since the first afterglow in 1997, the num¬ 
ber of detected hosts with spectroscopic redshift is only 


1/4 (374) of the total number of identified afterglow^ 
The majority are long GRBs (351), 23 are shorl[^ Eigh¬ 
teen years after the first afterglow localization, we are still 
dealing with a small number statistics. Nevertheless, the 
fact that GRBs are distributed over the entire redshift in¬ 
terval ever explored in the history of human kind (Eig.[^, 
from ^ = 0.0085 (GRB 980425 at 37 Mpc from the Milky 
Way; Galama et al. 1998) to z = 8.23 (GRB 090423; Sal- 
vaterra et al. 2009; Tanvir et al. 2009a), makes them the 
most valuable resources of exploration of the dawn of the 
Universe. 

2. The typical GRB host galaxy 

The definition of the population of galaxies hosting 
GRBs is not precisely defined. When an optical after¬ 
glow is detected and a spectrum is obtained, it is often 
possible to measure the redshift from the identification of 
absorption lines, for instance the strong MgHAA2796, 2803 
doublet (in the optical for the wide redshift interval 0.35 < 
z < 2.5). We automatically assume that these absorption 
features are associated with the interstellar medium (ISM) 
of the host galaxy, thus implicitly conclude that the host is 
identified. However, when we talk about the properties of 
the host galaxy population as a whole, we generally refer 
to the detected stellar and gas emission of the galaxy which 


^However, we notice that the fraction of detected hosts is much 
higher and at least 80% for well defined complete samples (e.g., 
TOUGH, or SHOALS surveys). 

^This large long-to-short number ratio is partly and likely due to 
the rarity of short events, but mainly determined by the faintness of 
their X-ray/optical afterglow, which, together with the short dura¬ 
tion, makes a precise localization and identification more difficult 
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Figure 2: Comparison between two high- 2 ; GRB afterglow spectra 
(probing the cold ISM in the host galaxy). Upper panel: spectrum 
of the afterglow of GRB 090323 at 2 ; = 3.57 (Savaglio et al. 2012), 
with strong metal absorbers revealing the presence two systems sep¬ 
arated by 660 km s“^ (total A^hi = with super solar 

metallicity ([Zn/H] = +0.25 ± 0.09). Lower panel: spectrum of the 
afterglow of GRB 090926 at 2 ; = 2.1062 (Rau et al. 2010), character¬ 
ized by a strong HI absorption (A^hi = weak metal 

lines, from which a low metallicity is derived ([Si/H] = —1.89 + 0.17). 
At 2 ; = 3.57, the Universe was 1.8 Gyr old and at 2 ; = 2.11, it aged 
by 1.36 Gyr. 

lies closest to the afterglow. By doing this, we not always 
give the proper emphasis to the results obtained from ab¬ 
sorption line studies, because the direct identification of 
the galaxy is not always possible. Our view becomes more 
incomplete at ^ > 1.5, where galaxies become increasingly 
fainter. As a consequence, the ’GRB host population’ 
we generally discuss includes basically galaxies detected 
in emission, photometrically and spectroscopically, mainly 
at 2 ; < 1.5, spanning about 2/3 of the entire history of the 
Universe. What about the GRB hosts during the first 4.3 
Gyr of life? 

For only less than half of all GRBs with measured red- 
shift (160) imaging has shown the presence of a galaxy. 
Multi-band photometry from optical to the near infrared 
(NIR) at low redshift is sufficient to derive with good 
accuracy the stellar mass of the galaxy. Of this popu¬ 
lation, another half or so are those with one (or more) 
emission lines detected, typically the strong [OIIJA3727 
feature, more easily detected (in the optical/NIR) when 
z < 1. Emission lines provide a direct measurement of 
SFR, the most known parameter. If, together with the 
oxygen, one Balmer line is detected, then the metallicity 
can be estimated. This can be H/d for ^ < 1, otherwise 
Ra for z < 0.5. At the moment, for at most 30 galaxies it 
was possible to determine the metallicity. When Ra and 


R/3 are detected together, the line ratio (Balmer decre¬ 
ment) gives the dust extinction in the optical Ay (about 
22 hosts). These samples are small, but we can improve 
the statistics by considering gas properties from absorp¬ 
tion lines. 

At high redshift, GRB afterglow spectroscopy provides 
independent measurements of the metallicity (at least 26 
GRB hosts). This is possible exclusively for 2 ; > 2, when 
the Lyman-a absorption line (necessary to measure the 
HI column density) is redshifted to the optical window. 
This redshift limit is set by the fact that the only suit¬ 
able UV spectrograph, on Hubble Space Telescope, is not 
used to observe GRB afterglows. If HI is not detected, 
low-ionization lines, associated with heavy elements in the 
neutral ISM, are used to derive relative abundances. For 
instance, the zinc-to-iron or silicon-to-iron ratio are di¬ 
rectly connected to the amount of dust depletion, there¬ 
fore, to the dust column distributed along the sight line in 
the host. 

The overlap between these different samples is still 
very limited. For only a very few individual GRBs it 
was possible to detect the afterglow and study absorption 
lines, together with the host in imaging, and also perform 
the spectroscopy for the emission lines, therefore, derive 
in principle the maximum number of galaxy parameters. 
However, soon the heavy use of the Optical-NIR spectro¬ 
graph X-Shooter (at the Very Large Telescope) for GRB 
programs will deliver first results on relatively large num¬ 
ber^ which will allow the exploration of an unknown ter¬ 
ritory, in a way that even normal galaxy surveys cannot 
do. We have to emphasize that the scientific community 
dealing with GRBs cannot really be accused of being lazy 
for not publishing data and results. An interesting statis¬ 
tics is the number of papers where the information on GRB 
host galaxies is reportecQ more or less 450 are the papers 
from which information for about 250 GRB hosts are re¬ 
ported. Yet, our understanding of this galaxy population 
is effected by the small number statistics. The difference 
with normal galaxy surveys, that count many thousands of 
galaxies, is apparent. For instance, the Sloan Digital Sky 
Survey (SDSS), 15 years after the beginning of the project 
includes photometry and spectroscopy for a sample of half 
a billion galaxies. 

For a long time and still today, the scientific community 
has adopted the idea that the typical GRB host galaxy 
is a small star-forming metal poor galaxy with low dust 
content (e.g., Levesque et al. 2010a; Graham & Fruchter 
2013; Vergani et al. 2014). For instance, a typical host is 
the one of GRB 011121 (Bloom et al. 2002; Kiipcu Yoldas 
et al. 2007), among the first ever studied in detail. The 
galaxy is at low-redshift {z = 0.362), small (stellar mass 


^ At the time of writing, the results of a large sample of X-Shooter 
spectra of GRB hosts at 0.1 < 2 : < 3.6 are in the process of being 
published (Kriihler et al. 2015) 

"^See the GHostS database at http://www.grbhosts.org for more 
information 
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Figure 3: Fraction of column density of molecular hydrogen = 2N(H2)/[N(H1) + 2N(H2)] vs. the total column density of hydrogen 
A/’(HI) + 2A/’(H2). The blues dots, red triangles and green diamonds correspond to H 2 detections in the Small Magellanic Cloud (SMC, 
Tumlinson et al. 2002), Large Magellanic Cloud (LMC, Tumlinson et al. 2002), and Milky Way (MW, Savage et al. 1977), respectively. Black 
filled and empty squares are measurements and upper limits in QSO-DLAs, respectively (Noterdaeme et al. 2008; Cuimaraes et al. 2012; 
Muzahid et al. 2015). Pink filled stars are measurements and upper limits in GRB-DLAs (Fynbo et al. 2006; Tumlinson et al. 2007; Prochaska 
et al. 2009; Kriihler et al. 2013; Friis et al. 2014). The region below the dashed-red line indicates a column density of molecular hydrogen 
below logA/’(H 2 ) = 14.7. 


Mq), with low metallicity log{Z/ZQ) = 
— 1.16 and low dust extinction {Ay = 0.4), but, given its 
mass, a relatively high SFR^ 2 Mq yr“^. The initial view, 
according to which the host of a GRB, whose progenitor 
is a massive and metal poor star, is special and does not 
represent the general galaxy population, has been hard to 
change. 

These common properties, derived at low redshift, can 
naturally be explained if we consider that galaxies are eas¬ 
ier to observe and study if they are close, at ^ < 1, where 
faint targets can be detected, emission lines measured in 
the optical, and dust extinction, SFR and metallicity esti¬ 
mated. At this low redshift, the most common star form¬ 
ing galaxy is small. High star formation activity is the 
necessary ingredient to have a high chance of a massive- 
star explosion. Small galaxies are typically metal and dust 
poor. The natural consequence is that the host galaxy of a 
low-z GRB is a small galaxy which is also metal and dust 
poor. However, extrapolating this idea to any redshift is 
itself a bias. Instead, it should be more natural to think 
that GRBs at high redshift happen in galaxies different 
from does at low redshift, because the Universe back then 
was very different. 

It would take too long to describe the large variety 
of properties seen lately in GRB host galaxies. We re¬ 
port a few examples. For instance, the two cases seen in 
Fig-i where the two afterglow spectra for GRB 090926A 
at ^ = 2.11 (Ran et al. 2010) and GRB 090323 at ^ = 3.57 
(Savaglio et al. 2012) are compared. Their metallicity are 
at two extremes: the former is about 1/100 the solar value 
(large HI absorption and the weak metal lines), the latter is 
over solar (strong metal lines and relatively low HI absorp¬ 


tion). It is remarkable in the latter that a high chemical 
enrichment is measured in a galaxy when the Universe was 
less than 2 Gyr old. Unfortunately the host is detected, 
but its mass is not known. 

3. The meaning of the observational bias 

Since the first GRB redshift measurement in 1997, many 
scientists inside and out the GRB community have learned 
that galaxies hosting GRBs do not represent the whole 
galaxy population. It is believed that GRBs are special 
and rare events and their hosts are observational biased. 
Moreover, theoretical models predict that, for a GRB to 
occur, a jet of relativistic particles has to break from the 
star collapsing core. The gamma-ray radiation has to 
emerge from inside and stay collimated, because if the 
emission were isotropic, it would be too large (equivalent 
to the rest-mass energy of the star). The collimation is pre¬ 
served by the high angular momentum, which is possible if 
radiation pressure (less efficient in a low-metallicity star) 
is low, and the mass loss as well. Therefore, it is easier to 
imagine that the GRB ambient is chemically pristine. 

However, it is strange that the dozen spectroscopically 
identified SNe, seen within 2 weeks after the GRB (mainly 
^ < 0.5), are type Ic: no hydrogen nor helium, but silicon. 
These features are explained by mass loss of the external 
envelope, through stellar wind (like in Wolfe-Rayet stars) 
or by the presence of a companion in a binary system. 

From one side, this issue is fundamental and still highly 
debated. From the other, it was clear since the beginning 
of the first redshift that there was something special about 
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Figure 4: Dust extinction in the visual band Ay derived from GRB afterglows and QSO-DLAs. Blue circles are Ay derived from the SED 
fitting of GRB afterglows from X-ray to UV and optical (Greiner et al. 2011; Schady et al2011; Zafar et al. 2012). Filled triangles are the 
same, but using UV to optical only (Kann et al. 2006; 2010; 2011). The green histogram indicates the mean Ay and dispersion per redshift 
bin from a large sample of QSO-DLA, derived by using the zinc-to-iron relative abundances and assuming that these are proportional to 
the visual dust extinction Ay (see Savaglio 2006). The black curve indicates the extinction in the rest-frame U-band for a constant U-band 
extinction in the observed frame Ay =2.7. 


the host galaxies. They looked different, although not nec¬ 
essarily unique, from the observational point of view. It 
looked like we had missed these galaxies while using the 
traditional approach: observational campaigns start from 
optically bright targets. With GRBs, we can in principle 
overcome this limitation: GRBs are bright in the gamma- 
ray and short, and pinpoint to locations in galaxies re¬ 
gardless of their brightness. It is hard to imagine a less 
biased way of selecting galaxies. The main problem of this 
population right now is, as said above, the small num¬ 
ber statistics, and not the way (still vaguely understood) 
GRBs select galaxies. 

The GRB host population suffers of some kind of bias, 
but not more than any other galaxy population we can 
think of. Any galaxy population is detected with specific 
tools and using some selection criteria, and by definition 
any survey is biased. Observational biases are corrected by 
using well tested and complex analysis methods, but this 
is done by assuming that we have full control of the ob¬ 
servational tools (possible) and have an idea (up front) of 
what we are missing. GRB hosts are selected in a very dif¬ 
ferent way, and this, by itself, can only be good because it 
tackles the problem of galaxy search from a different point 
of view. Being different is not a good reason for discarding 
them when studying galaxy formation and evolution. This 
would introduce a bias, and there is no scientific justifica¬ 
tion to do so. 

GRB host galaxies have an important value for our 
understanding of galaxies at any epoch. In the 450 pa¬ 
pers mentioned above, a large variety of telescopes, wave¬ 


length bands, instruments, depths and analysis methods 
are used. Thus, very different are the properties that we 
found, which is by itself an indication that, given the small 
number statistics, there cannot be a typical GRB host 
galaxy. The most common comments the GRB commu¬ 
nity faces when applying for telescope time is that ”the 
sample is too small” or ’’biased and does not represent 
the galaxy population”. This is not always a good justi¬ 
fication: we have achieved many fundamental discoveries 
from a small sample, adding one or a few more objects is 
important in any case. 

4. Spanning the entire history of the Universe 

There cannot be a typical GRB host simply because 
GRBs are seen from z = 0 to the highest ever measured 
redshift ^ = 8.2, when the first collapsed objects formed. 
Over this time, the Universe has dramatically changed. 

Not surprising if GRB hosts did as well. It is widely ac¬ 
cepted that the Universe was much more active in the 
past, as clearly seen from the redshift evolution of the star- 
formation rate density (SFRD). Fig.shows the history of 
the SFRD, characterized by a beginning, a steady/intermediate 
state, and a decaying era. Our knowledge is much more ac¬ 
curate for the last 8 Gyr {z < 1), when the SFRD dropped 
by at least a factor of 10. This is the Universe where half 
of all GRB hosts are detected and from where most of our 
knowledge comes from. 

At higher redshift, hosts are harder to see in emission 
because they are fainter, but the GRB can overcome this 
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Figure 5: Comparison of two composite UV spectra probing the cold ISM in galaxies. The red spectrum is obtained from ~ 25 GRB afterglow 
spectra in the redshift interval z = 0.9 — 1.5 (Christensen et al. 2011). The black one includes about 100 QSO-DLA spectra in the redshift 
range 2 < z < 4 and is characterized by large HI column densities Nm > 5 x 10^^‘° cm“^ and relatively low metallicity ZjZ^ ~ 1/20 
(Noterdaeme et al. 2014). Vertical marks indicate resonance (solid blue) and fluorescent (dashed red) lines. The GRB-DLA composite is 
characterized by much stronger metal lines than the QSO-DLA one. 


limitation. Afterglow spectra reveal absorption lines as¬ 
sociated with the ISM in the host. Column densities of 
the gas are often large, and the absorbers are classified 
as damped Lyman-a systems (DLAs). In these cases the 
gas is mainly neutral, the ionization correction can be ne¬ 
glected and the metallicity can easily be derived. This 
was done for GRBs in the interval 2 < z < 6.3. At z > 4, 
the redshift is measured for 27 GRBs from the afterglow 
(Fig. 0 , but not much luck for the direct detection of the 
host: only 4 have been seen, in the rest-frame UV, and 
none at 2 : > 5 (where 10 are the total number of GRB 
redshifts). 

GRB-DLAs have often a large column density of HI. Of 
the 62 GRB-DLAs with a measured Vhi > 2 x 10^^ cm“^, 
53 (85.5%) have Nm > and 13 (21%) have 

A^hi > 10^^ cm ^. The largest value ever measured {Nm = 
1022-70±o. 15 (3m“2) was found in GRB 080607 dit z = 3.0363 
(Prochaska et al. 2009). For comparison, the fraction of 
DLAs with Nm > cm“2 found along QSO sight lines 
is 5.4% only. This is one of the most striking differences 
between what has been known for decades about the ISM 
in high- 2 ( galaxies and what is now observed with GRB 
spectroscopy. 

Another science topic where GRB investigations re¬ 
vealed interesting results is the molecular gas. The low 
fraction of molecular hydrogen in QSO-DLAs is a well 


known phenomenon. Fig. 3 shows the H 2 column densities 
measured in local environments (the Milky Way and Mag¬ 
ellanic Clouds, Savage et al. 1977; Tumlinson et al. 2002), 
and in QSO-DLAs. Many are the low values and up¬ 
per limits in the former population. The investigation of 
molecular hydrogen was possible only though a very lim¬ 
ited number of GRB afterglows, eight in total, and half 
of them led to a positive detection. Remarkable is the 
relatively high fraction of molecular hydrogen fog/H 2 = 
2V(H2)/[V(HI) + 2V(H2)] - -1.14 found in GRB 080607 
at 2 : = 3.0363 with the extreme column density of neutral 
hydrogen Nm = 10^^-^ (Prochaska et al. 2009). 

With a better statistics, available in the last few years, 
many are the GRB hosts revealing unexpected (according 
to the common belief) properties. A few have high molec¬ 
ular gas content (see above), others are massive, many 
have a large amount of dust and dust extinction. The 
dust extinction is easily estimated by assuming that the 
intrinsic emission is a simple or broken power law, and fit¬ 
ting the observed SED of the afterglow in NIR, optical, 
UV and X-ray. In Fig. 4, the optical extinctions Ay in 
GRB afterglows is shown together with the estimated av¬ 
erage (plus dispersion) Ay in QSO-DLAs. In the latter. 
Ay is estimated by considering that the dust depletion, 
represented by the zinc-to-iron relative abundance, gives 
the dust content along the sight line, thus, the dust extinc- 
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Figure 6: Examples of host galaxies associated with dark GRBs, and detected with Herschel (Hunt et al. 2014). Filled green circles are the 
first ALMA detections of GRB hosts (Hatsukade et al. 2014). 


tion (Savaglio 2006). The figure shows that GRBs probe 
dustier environments than QSO-DLAs. 

The presence of dust is also manifested by the UV 
bump, the broad absorption feature centered at A ^ 2175 
A. Prominent in the Milky Way, it is much weaker in the 
LMC and absent in the SMC. Its origin is still not to¬ 
tally known: graphites and silicates, mixed with complex 
molecules? The UV bump is not easily detected outside lo¬ 
cal environments. In QSO-DLAs, it is generally very weak 
and can be seen only from composite spectra (see Ledoux 
et al2015 for a more detailed discussion). 

A strong UV bump is seen in several GRB afterglow 
spectra. For instance, Prochaska et al (2009) detected it 
in GRB 080607 at 2 : = 3.036, from which Ay = 3.2 was 
derived. In the same spectrum, a large HI column density 
and molecular absorption is measured as well (see Section 
[^. The host of this dark GRB was widely investigated 
from the UV to the infrared, and a relatively high stellar 
mass is derived, 4 x 10^^ Mq (Perley et al. 2013). 

This is another host which, according to the canonical 
view, can be considered ’abnormal’. More recently, an 
unusually strong UV bump is identified in GRB 140506 
at 2 : = 0.889. The best-fit model gives an extreme Ay 
(Fynbo et al. 2014). In general, the presence of the UV 
bump is also demonstrated from multi-band photometry 
and SED fitting of the afterglow. One clear detection is in 
GRB 070802 at z = 2.45, giving Ay = 1.8 and a best-fit 
LMG extinction law (Kriihler et al. 2008). 

5. The properties of GRB host galaxies 

Not all GRB hosts are metal poor. At ^ < 2, a few 
exceptions were found (Levesque et al. 2010b; Perley et al. 
2012; Kriihler et al. 2012; Niino et al. 2012). At 2 : > 2, 
GRB-DLAs display a large dispersion. We already men¬ 
tioned the two extreme measurements in Fig.|^ The super¬ 
solar metallicity found in GRB 090323 at 2 : = 3.57 (Savaglio 


et al. 2012) demonstrates that even at high redshift, the 
host is not necessarily metal poor. 

The typical high metallicity in GRB-DLAs is more 
clearly seen in Fig. 5. This displays an instructive compos¬ 
ite spectrum in the rest-frame interval AA = 2200 — 2900 A 
(around the Fell and Mgll absorptions) from ^ 25 GRB 
afterglow spectra (Ghristensen et al. 2011). The redshift 
interval is from 2 : = 0.9 (blue end of the spectrum) to 
^ = 1.5 (red end of the spectrum). The metallicity is not 
measured in the sample, but cold-ISM absorption lines are 
strong, especially in optically dark bursts. Dark bursts are 
mostly not seen in the optical and have an optical to X- 
ray spectral index /^ox <0.5. The global full composite 
includes the spectra of 60 GRB afterglows, mainly in the 
redshift interval 2 < ^ < 2.4, with an average metallicity 
which is relatively high: ZjZ^ ^1/6. In the same fig¬ 
ure, we see the comparison with the composite spectrum 
of gas-rich QSO-DLAs at 2 < 2 : < 4, characterized by 
Vhi > 5 X 10^^’^ cm“^ and relatively low metallicity, with 
ZjZ^ ^ 1/20 (Noterdaeme et al. 2014). 

A detailed comparison with more composite spectra, 
representing a large variety of galaxy populations, is given 
in Savaglio (2013). The average spectrum of 13 massive 
galaxies at ^ ^ 1.6 (median stellar mass = 2.4 x 10^^ 
Mq, SFR = 30 Mq yr“^, specific star-formation rate sSFR 
= 1.2 Gyr“^; Savaglio et al. 2004) shows similar absorp¬ 
tions to the GRB composite in Fig. 5. The composite 
of a complete sample of UV-bright z ^ 1 galaxies has 
much weaker absorption lines, but sizable emission lines 
(from the hot gas), with the tendency of stronger absorbers 
to be more common in brighter galaxies (Martin et al. 
2012). These galaxies have SFR = I — 100 MQyr“^ and 
M* = 10^-^ - lO^i-^ MQyr-i {sSFR = 0.07 - 6 Gyr-^). 
A similar composite is the one of 28 local {z < 0.05) 
starburst and star-forming galaxies, with median metal¬ 
licity \ogZ/ZQ = —0.5, UV luminosity and LC-band abso¬ 
lute magnitude I /1500 = 5 x 10^^ ergs“^ and Mk = 
—21.35, respectively (Leitherer et al. 2011). Using the 
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Figure 7: Observations of the field of GRB 090423 (z = 8.23). Left, center and right: HST (Tanvir et al. 2012), Spitzer (Berger et a. 2014) 
and ALMA (Berger et al. 2014) images. The scale is the same for all images. Ten arcsecs ai z = 8.23 corresponds to a physical size of 48 
kpc. The position of the GRB (indicated by the cross and obtained from NIR imaging of the afterglow; Tanvir et al. 2009b) is known with a 
precision of 0.3 arcsecs (1.45 kpc). 


empirical relations in Savaglio et al. (2009), these num¬ 
bers are translated into a star-formation rate SFR1500 ^ 1 
M0yr“^ and a stellar mass ^ 6 x 10^ Mq (assuming 
Ay = 1). This gives a (rather uncertain) sSFR of a few 
Gyr“^. Surprisingly, these values are not very dissimilar 
from those of the 2: ^ 0.75 GRB host sample (Savaglio et 
al. 2009), despite the apparent difference with the GRB 
composite. However, we notice that the median redshift 
of the GRB host sample is lower than the redshift interval 
covered by the GRB composite {z = 0.9 — 1.5), indicating 
again a redshift evolution of the galaxy population hosting 
GRBs. 

6. More massive GRB hosts at 2 : > 1.5 

As described above, at 2; > 1.5, some GRB hosts are 
metal rich, massive, dusty (dark GRBs), or highly star 
forming (Hunt et al. 2011; Kriihler et al. 2011; Rossi et al. 
2012; Hunt et al. 2014). This kind of galaxies can be very 
bright in the sub-millimeter. A few examples, detected 
with Herschel and ALMA, are shown in Fig. 6 (Hunt et 
al.2014; Hatsukade et al2014). However, a systematic 
search in a large sample at 2: < 1 shows that the total 
number of radio bright hosts is very small (Michalowski 
et al. 2012; Perley & Perley 2013). One possible explana¬ 
tion is that sub-mm galaxies (SMGs) account for at most 
20% of the cosmic SFRD (Michalowski et al. 2010). Fu¬ 
ture surveys can bridge the gap at 2 < ^ < 4 and explain 
a possible steep redshift evolution. 

The fraction of pair absorbers in ^ > 1.5 GRB af¬ 
terglow spectra has been found to be almost three times 
higher than in QSO-DLAs (which probe random galaxies), 
suggesting that galaxy interactions may play a role in the 
formation of massive stars at high redshift (Savaglio et al. 
2012). Another indication is the large fraction (at least 
40%) of known GRB hosts at 2; > 1.5 showing interac¬ 
tion, disturbed morphologies, or galaxy pairs (Thone et 


al. 2011; Vergani et al. 2011; Chen 2012; Kriihler et al. 
2012). The interaction hypothesis is not surprising if one 
considers the higher fraction of galaxy mergers seen in the 
past of the Universe with respect to today (Bluck et al. 
2012). 

7. High-redshift GRB hosts 

The investigation of GRB hosts at 2: > 4 has been 
so far particularly difficult. At 2; > 4 and 2; > 5, we 
know 27 and 10 GRBs with measured redshift, respec¬ 
tively. Intensive search of 4 < 2; < 5 GRB fields resulted 
in rest-frame UV (e.g., star formation) detections for 7 
hosts. Detailed analysis of the luminosity function for de¬ 
tected (and non-detected) hosts at 3 < 2; < 5 suggests that 
these galaxies might represent the most known population 
of galaxies at those redshifts, the Lyman break galaxies 
(Greiner et al. 2015). However, at 2: > 5, 5 GRB fields 
have been observed (Basa et al. 2012; Tanvir et al. 2012, 
Berger et al. 2014), with negative results. The lack of Ly¬ 
man break galaxies at very high redshift is also found by 
Schulze et al. (2015), though results are not totally consis¬ 
tent with those in Greiner et al. (2015. If no dust correction 
is assumed (dust is not expected to be abundant in a less 
than 1.2 Gyr old Universe), the UV-luminosity limit I/1500 
can be translated into a SFR limit (Savaglio et al. 2009). 
The deepest search obtained with HST/NIR (Tanvir et 
al. 2012), Spitzer medium IR, and ALMA radio observa¬ 
tions (Berger et al.2014) for the host of GRB090423 at 
^ = 8.23 gave negative results (Fig. 7). The limiting mag¬ 
nitude obtained with HST in the NIR is itiab > 30.29, 
which corresponds to a very low upper limit for the UV 
luminosity, I/1500 < 3.7 x 10^^ erg s“^, and an exception¬ 
ally low SFR < 0.07 M0yr“^. 

We can compare these low SFR limits to the stellar 
mass expected from numerical simulations. About 70% of 
the hosts at 2; > 6 predicted by Salvaterra et al. (2013) 
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have stellar mass in the range = 10^ — 10^ Mq, while 
the star-formation rate and metallicity are in the intervals 
SFR = 0.03 — 0.3 Mq yr“^ and ZjZ^ = 0.01 — 0.1, respec¬ 
tively. The SFR limit for the host of GRB 090423 indicates 
that a very low stellar mass, ^ 10^ Mq, is possible if 
sSFR < 70 Gyr“^. Vice versa, if we assume sSFR ^ 10 
Gyr“^, the limit SFR < 0.07 Mq yr“^ gives < 7 x 10^ 
Mq. We conclude that in the past, GRB hosts must have 
been very small. 

8. Are some high- 2 ; GRBs hosted by young star 

clusters? 

How about if a fraction of GRBs at high redshift hap¬ 
pened in young progenitors of today’s globular clusters? 
If so, finding a host like that might still be out of present 
telescope capabilities. The upper panel of Fig.[^ shows 
a beautiful HST image of a well known globular cluster 
(GC): Messier 80 (M80). Due to its intrinsic brightness 
(absolute magnitude ^ — 7) and proximity to the Sun (10 
kpc), this is one of the most studied GCs ever. M80 is 
characterized by the usual properties for GCs: old (12.5 
Gyr), metal poor (log ZjZ^ —1-5), compact, with a high 
stellar density, and is also relatively massive (M 5 x 10^ 
Mq). If we believe that the Universe today is 13.7 Gyr old, 
M80 must have formed when the Universe was 1.2 Gyr old. 
That means that GRBs at z 5 would be contemporary 
to the newly born M80. However, one main problem of a 
proto-GC as a possible host of high-z GRBs is the small 
fraction of stellar mass contained in GCs with respect to 
the total stellar mass in galaxies. Less than 200 are the 
known GCs in the MW. If we assume they all have a mass 
similar to M80 (which is a massive GC) and that the stellar 
mass of the MW is M^{MW) 4 x 10^^ Mq, the fraction 
of stellar mass in GCs with respect to the total cannot be 
higher than M>^(GCs)/M>^(MW) ~ 0.25%, very low. 

However, the picture cannot be discarded too quickly. 
Young proto-GCs must have had very different proper¬ 
ties, perhaps similar to a young massive star cluster to¬ 
day. Moreover, it has been suggested that intense star 
formation in starbursts create the physical conditions to 
form numerous and massive star clusters (e.g., Adamo et 
al. 2011). Finally, hydro dynamical simulations have shown 
that galaxy collision has a role in the formation of massive 
star clusters (Renaud et al. 2015), which eventually will 
evolve in GCs. As mentioned in Section 6, interactions 
and mergers are common in high-2; GRBs. 

The lower panel of Fig.shows the recently discovered 
massive star cluster R136, in the star-forming region 30 
Doradus of the LMC (Crowthler et al. 2010). This star 
cluster is young (a few Myr old), and as massive as M80. 
What is interesting from the point of view of GRBs is that 
R136 contains a high concentration of massive stars, over 
70 O and Wolfe-Rayet stars within 5 pc. Equally impor¬ 
tant, Conroy (2012) found indications that the initial mass 
of GCs were at least 10 — 20 times larger than today, eas¬ 
ily exceeding 10^ Mq. Mass loss is due to stellar winds. 



Figure 8: Top: HST/WFPC2 image of one of the densest globular 
clusters (GCs) known in the Milky Way: Messier 80 (credit: The 
Hubble Heritage Team - AURA/ STScI/ NASA). The estimated 
mass, metallicity, and age are M ~ 5 x 10^ Mq, \og Z/Zq = —1.47, 
and t ~ 12.5 Gyr, respectively. Its age corresponds to a redshift of 
formation 2 : ~ 5. At that time, this GC must have been dominated 
by young and massive stars, with many of them being in binary or 
multiple systems. Bottom: HST UV, V and R-band image of the su¬ 
per star cluster R136, in the region 30 Doradus of the LMC (credit: 
R. O’Connell, University of Virginia, and the WFC3 Science Over¬ 
sight Committee). Similar massive star clusters are the progenitors 
of GCs, and might host a GFB during the present active phase. 
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shocks and SN events. Moreover, a certain number of GCs 
in the MW might have disappeared early on, being grav¬ 
itationally unbound within 1 Gyr (Vesperini et ah 2009), 
or disrupted during the merging of proto-galaxies. 

From the Galaxy point of view, it has been known for a 
decade now that a large fraction of the cosmic stellar mass 
was produced at 2; > 1, when the Universe was almost half 
of its present age. In a recent review article, Madau & 
Dickinson (2014) show that the stellar-mass density from 
z = b to now has increased by a factor of ~ 60 (Fig. 9). A 
galaxy like the MW was still in the process of forming at 
^ = 5, very likely the main body of its mass (for instance, 
the disk) was not in place jet, whereas likely a proto-bulge 
was. However, a large fraction of the proto-GGs, those 
that today are older than 12.5 Gyr were already formed 
or on the process of forming, some will never survive for 
long, suffering from the disruptive gravitational force of 
the MW. 

All together, this means that the initial fraction of stars 
in proto-GGs must have been much larger than the 0.25% 
fraction of today. We quantify under which circumstances, 
at z = 5, the total mass in proto-GGs would be about the 
mass of the proto-MW. If we assume that: (a) proto-GGs 
were 10 times more massive, which means M(proto-GG) ^ 
5 X 10^ Mq (6) the total number of proto-GGs was lower 
than today because many formed at 2: < 5; (c) but this 
is compensated by assuming that many first generation 
proto-GGs did not survive until today; (d) the stellar mass 
of the MW was 40 times lower, e.g., M(proto-MW) ^ 
1 X 10^^ Mq, then all these conditions together give a ratio 
of stellar mass in proto-GGs to the one of the proto-MW 
M*(proto-GGs)/M*(proto-MW) ^ 1. Stars in proto-GGs 
and the proto-MW are equally important. 

Another last situation in favor of proto-GGs as hosts 
of high-^ GRBs is that they have a high density of bi¬ 
nary systems. Massive and rapidly rotating binary sys¬ 
tems are expected to be very common in a low-metallicity 
environment (Sana et ah, 2012) and such a configuration 
is considered to be very favorable for high-redshift GRBs 
(Bromm & Loeb(2006). For the description of the forma¬ 
tion of fast-rotating massive stars in a possible scenario of 
proto-GGs, see Krause et al. (2013). 

In Fig. 10, we show this possible scenario, with a GRB 
happening in a proto-GG at 2: = 5. On top an optical 
image of the Andromeda galaxy. In the center a UV image 
of the same galaxy, with artificially overlaid many GGs like 
M80. At the bottom, we imagine that at 2: = 5, only the 
bulge of such a giant spiral galaxy was in place, the UV 
part in the disk with young stars today would form later, 
while all GGs were already there. In the same image, a 
GRB has exploded inside a GG. 

9. Conclusions 

The impact of GRB host galaxies on the understand¬ 
ing of galaxy formation and evolution is still affected by 


lookback time (Gyr) 


0 2 4 6 8 10 12 



Figure 9: The evolution of the stellar mass density in the Universe 
(from Madau & Dickinson 2014). The data points are collected from 
papers published since 2007. The solid line shows the prediction gen¬ 
erated by integrating the instantaneous star formation rate density 
as a function of redshift. The increase in the cosmic stellar mass from 
2 ; = 5 to today was a factor of ~ 60. Most of the action took place 
at 2 ; > 1. Below this redshift (shaded area), the Universe (which was 
less than 6 Gyr old, or 43% of its present age) had already formed 
more than 60% of the stars today. 

small number statistics. Their knowledge is mainly lim¬ 
ited to the 2: < 1.5 regime, where about half of the galax¬ 
ies hosting GRBs have been found and studied in detail. 
From this, the most accepted picture is that GRB hosts 
are generally small, star forming and metal and dust poor 
objects. However, at ^ > 1.5, metallicity, mass and dust 
extinction show a large spread, suggesting a different pop¬ 
ulation. Many host galaxies show disturbed morphologies, 
interactions with nearby galaxies and mergers. All this 
is nicely connected to the idea that local massive ellipti¬ 
cals today were young and bursty in the past, with some 
of them experiencing close encounters with other galaxies, 
which likely triggered intense episodes of star formation. 

At very high redshift, 2: > 5, the situation might have 
changed again. Massive galaxies were very rare, but these 
are observationally easier to find. In fact, at these dis¬ 
tances, deep searches failed to detected any GRB host, 
and relatively low SFRs were inferred. Unless dust con¬ 
tent was already important back then (unlikely), low SFRs 
means low galaxy mass. Therefore, GRB hosts in the past 
could have been more similar to the local counterparts. 
This is supported by recent findings according to which 
the UV fiux that ionized the Universe at high redshift was 
produced by very small, young and star forming galaxies 
(Bouwens et al. 2015). 

As an alternative picture, we propose a different paradigm, 
never considered in detail, according to which a proto 
globular cluster, which would look like a massive young 
star cluster in the local Universe, could host high-2: GRBs. 
Globular clusters today contain a very low fraction of the 
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UV image of Andromeda vvith M8p’s 


GRB 




Andromeda at z = 5 and ^ GRB in primordial M80 



Figure 10: A possible scenario for the origin of 2 ; > 5 GRBs. Top: 
optical image of Andromeda; at center: UV image of the same galaxy, 
with artificially overlaid several globular clusters as M80; bottom: the 
imaginary proto-bulge of the same galaxy, no disk present yet, and 
the same much younger proto-GCs. In one of the them, the explosion 
of a massive star in a binary system became a GRB. 


total stellar mass, but at z > 5 proto-GCs can be equally 
important in terms of stellar mass because many lost a 
large fraction of the mass in subsequent encounters with 
more massive galaxies. Moreover, giant galaxies like the 
MW have formed a large fraction of their present stellar 
mass in recent times. In support of this idea is the fact 
that proto-GCs probably contains a large fraction of mas¬ 
sive and rapidly rotating binary systems, also a favorable 
situation for GRBs. 

In summary, GRB hosts z < l.b are generally small 
galaxies, whereas at intermediate redshift an important 
fraction of hosts are massive, dusty or metal rich galax¬ 
ies. Going further in the past, in the primordial Universe, 
host galaxies are again likely very small systems, so small 
that, perhaps, these systems are massive young star clus¬ 
ters with a large fraction of massive and rotating binary 
systems, possible seeds of GRBs. 
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